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Abstract 

We develop a comprehensive quantitative description of the cross-section mech- 
anism discovered several years ago by Lazarian. This is one of the processes that 
determine grain orientation in clouds of suprathermal cosmic dust. The cross- 
section mechanism manifests itself when an ensemble of suprathermal param- 
agnetic granules is placed in a magnetic field and is subject to ultrasonic gas 
bombardment. The mechanism yields dust alignment whose efficiency depends 
upon two factors: the geometric shape of the granules, and the angle $ between 
the magnetic line and the gas fiow. We calculate the quantitative measure of 
this alignment, and study its dependence upon the said factors. It turns out 
that, irrelevant of the grain shape, the action of a flux does not lead to alignment 
if <I) = arccos(l/\/3). 
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I. INTRODUCTION. THE PHYSICAL NATURE OF EFFECT. 



Polarisation of the starlight is a long- known effect. Due to the correlation of polarisation 
with reddening, this phenomenon is put down to the orientation of particles in the cosmic-dust 
nebulae (Hall 1949, Hiltner 1949). This orientation causes differential extinction of electromag- 
netic waves of different polarisations (effect known in optics as linear dichroism) , and provides 
a remarkable example of order emerging in a seemingly chaotic system. 

In a nutshell, the polarization is due to the fact that the grains are non-spherical (i.e., have 
different cross sections in the body frame) and these non-spherical cross sections are somehow 
aligned within the cloud. A remarkable feature of this phenomenon is that, whatever orien- 
tational mechanisms show themselves in the dust-grain rotational dynamics, the orientation 
always takes place relative to the interstellar magnetic field. 

Thus, whenever the word "orientation" is used, it always is an euphemism for "orienta- 
tion with respect to the magnetic line" . Another semantic issue is the conventional difference 
between the meaning of words "orientation" and "ahgnment" . Historically, the word "orienta- 
tion" has been used to denote existence of one preferred direction in a physical setting. The 
analysis of cosmic-dust orientation has demonstrated that in all known cases this preferred 
direction is equivalent to its opposite. Simply speaking, were an instantaneous inversion of the 
magnetic field in the dust cloud possible, it would not alter the orientation of the dust ensam- 
ble. Because of such an invariance, the word "orientation" is often avoided and substituted 
by the word "alignment" which is thereby imparted with the desired broader meaning. (The 
adjective "orientational" , though, remains in use.) We shall abide by this verbal code. 

The rotational dynamics of an interstellar particle is determined by a whole bunch of accom- 
panying physical processes whose combination produces a variety of orientational mechanisms. 
Which of these come into play in a particular physical setting, depends upon the suprather- 
mality of the dust cloud. Suprathermal are, by definition, grains which spin so rapidly that 
their averaged (over the dust ensemble) rotational kinetic energy < Erot > much exceeds the 
(multiplied by the Bolzmann constant k) temperature Tgas of the surrounding environment. 
The suprathermality degree is then introduced as the following ratio: 
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P = < Erot > /f^Tgas . (1.1) 

Dust ensembles with (3 not much different from unity are called thermal or Brownian. Clouds 
with (3^1 are called suprathermal. In the observable Universe, values of f3 of order 10^ 
are not unusual. 

The leading reason for suprathermal rotation is formation of H2 molecules at the defects 
on the granule surface: over such a defect (called active site), two atoms of H couple to form 
a molecule, ejection whereof applies an uncompensated torque to the granule surface (Purcell 
1979). These, so-called spin-up torques keep emerging at each active site until the site gets 
"poisoned" through the everlasting accretion. After that, some other active site will dominate 
the spin dynamics of the grain, by its H2 "rocket" . This change of spin state will, with some 
probability, go through a short-term decrease, to thermal values, of the grain's angular velocity. 
Such breaks are called "cross-overs" or "flip-overs". 

Another pivotal circumstance to be mentioned here is the existing evidence of paramagnetic 
nature of a considerable portion of dust particles (Whittet 1992) which makes them subject to 
the Barnett effect. This phenomenon takes place in para- and ferromagnetics, due to inter- 
action between the spins of unpaired electrons and the macroscopic rotation of crystal lattice 
(Stoner 1934). This coupling has its origin in the angle- dependent terms in the dipole-dipole 
interaction of neighbouring spins. It spontaneously endows a rotating para- or ferromagnetic 
body with a magnetic moment parallel to the angular velocity^] (Lazarian & Roberge 1997). 
Purcell offered the following illustrative explanation of the effect. If a rotating body contains 
equal amount of spin-up and spin-down unpaired electrons, its magnetisation is nil. Its kinetic 
energy would decrease, with the total angular momentum remaining unaltered, if some share 
of the entire angular momentum could be transferred to the spins by turning some of the un- 
paired spins over (and, thus, by dissipating some energy). This potential possibility is brought 
to pass through the said coupling. 

An immediate outcome from granule magnetisation is the subsequent coupling of the mag- 



Another contribution to the magnetisation comes from the electric charge carried by the granule. 
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netic moment M with the interstellar magnetic field B : the magnetic moment precesses 
about the magnetic line. What is important, is that this precession goes at an intermediate 
rate. On the one hand, it is slower than the grain's spin about its instantaneous rotation 
axis. On the other hand, the precession period is much shorter than the typical time scale at 
which the relative-to-B alignment gets established^. The latter was proven by Dolginov & 
Mytrofanov (1976), for magnetisation resulting from the Barnett effect, and by Martin (1971), 
for magnetisation resulting from the grain's charge. 

If we disembody the core idea of the Barnett effect from its particular implementation, 
we shall see that it is of quite a general sort: a free rotator, though conserves its angular 
momentum, tends to minimise its kinatic energy through some dissipation mechanism(s). This 
fact, neglected in the Euler-Jacobi theory of unsupported top, makes their theory inapplicable 
at time scales comparable to the typical time of dissipation (Efroimsky 2002). The needed 
generalisation of the insupported-top dynamics constitutes a mathematically involved area of 
study (Efroimsky 2000), which provides ramifications upon wobbling asteroids and comets 
(Efroimsky 2001), rotating spacecraft, and even precessing neutron stars (Triimper et al 1986, 
Alpar & Ogelman 1987, Bisnovatyi-Kogan & Kahabka 1993, Stairs et al 2000). Fortunately, 
in the case of cosmic-dust physics, we need only some basics of this theory. A free rotator has 
its kinetic energy minimised (with its angular momentum being fixed) when the rotation axis 
coincides with the axis of major inertia. In this, so-called principal state, the major-inertia axis 
Z , the angular velocity , and the angular-momentum vector J are all aligned. In other, 
complex rotation states, both the maximal-inertia axis Z and angular-velocity f2 precess 
about the angular momentum J . This precession is also called "wobble" or "tumbling", in 
order to distinguish it from the precession of the magnetic moment M about the magnetic 
line. Similarly, the wobble relaxation (i.e., gradual alignment of axis Z and of f2 toward 

^ A more exact statement, needed in Section V below, is that the period of precession (about B ) of 
the magnetic moment M (and of Z ahgned therewith) is much shorter than the mean time between 
two sequent flip-overs of a spinning granule (Purcell 1979, Roberge et al 1993) 
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J ) should be distinguished from the granule alignment relative to the magnetic field B : the 
latter effect is the eventual target of our treatise, while the former is merely a trend in the 
tapestry. Still, the wobble relaxation is far more than a mere technicality: it is important 
to know if a typical time of the wobble relaxation is much less than the typical times of 
the external interactions (like, say, the period of precession of M about B ). In case 
the wobble relaxation is that swift, one may assume that the precession (about B ) of the 
magnetic moment M is the same as precession of the angular momentum J about B : 
indeed, in this case, both J and the major-inertia axis Z will be aligned with O which 
is parallel to M . This parallellism of all four vectors is often called not alignment but 
"coupling", to distinguish it from the alignment relative to B . This coupling is enforced by 
two different processes. One is an effect kin to that of Barnett: tumbling of O relative to 
conserved^ J (and, therefore, relative to an inertial observer) yields periodic remagnetisation 
of the material, which results in dissipation. The other effect is the anelastic dissipation: in 
a complex rotational state the points inside the body experience time- dependent accelaration 
that produces alternate stresses and strains. Anelastic phenomena entail inner friction (which 
may be understood also in terms of a time lag between the strain and stress). The contributions 
from the Barnett and anelastic effects to the coupling were compared by Purcell in his long- 
standing cornerstone work (Purcell 1979). Purcell came to an unexpected conclusion that 
the input from the Barnett effect much outweights that from anelasticity. Having set out 
the calculations, Purcell continued with phrase: "It may seem surprising that an effect as 
feeble, by most standards of measure, as the Barnett effect could so decisively dominate the 
grain dynamics". After that, Purcell tried to pile up some qualitative evidence, to buttress 
up the unusual result. Still, the afore quoted emotional passage reveals that, most probably, 
Purcell's tremendous physical intuition signalled him that something had been overlooked 
in his study. An accurate treatment (Lazarian & Efroimsky 1999) shows that the anelastic 
dissipation is several orders of magnitude more effective than presumed, and in many physical 



At time scales shorter than a precession cycle about B . 
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settings it dominates over the Barnett dissipation. The case of suprathermal dust is one such 
setting. Without going into redundant details, we would mention that combination of the two 
dissipation processes provides at least partial alignment of Z and f2 toward J in Brownian 
clouds, and it provides a perfect alignment in suprathermal ones. 

The presently known mechanisms of grain alignment can be classified into three cathegories: 
mechanical mechanisms, paramagnetic mechanisms, and via radiative torques. The latter 
mechanism was addressed in (Dolginov & Mytrofanov (1976); Lazarian 1995b; Draine and 
Weingartner 1996, 1997). It has not yet been well understood. The paramagnetic alignment is 
due to the Davis-Greenstein (1951) mechanism (initially suggested for Brownian dust particles), 
and due to the Purcell (1979) mechanism (which is a generalisation of the Davis-Greenstein 
mechanism, to the suprathermal case). The Davis-Greenstein and Purcell processes operate to 
bring the granule's rotation axis (which is, as explained above, fully or partially aligned with the 
granule's major- inertia axis), into parallelism with the magnetic line. This happens because 
precession of grain's spin axis about B entails material remagnetisationQ and, therefore, 
dissipation resulting in a slow removal of the rotation component orthogonal to B. The induced 
alternating magnetisation M will lag behind rotating B, giving birth to a nonzero torque equal, 
in the body frame, to M x B. It can be shown (Davis & Greenstein 1951) that this torque 
will entail steady decrease of the orthogonal-to-B component of the angular velocity^. The 
so-called^ mechanical alignment comprises the Gold (1952) mechanism, and those of Lazarian 

^It is assumed that the grain is either paramegnetic (Davis & Greenstein 1951) or ferromagnetic 
(Spitzer & Tukey 1951), (Jones & Spitzer 1967). The case of a diamagnetic granule has not had been 
addressed in the Uterature so far. 

rigorous analysis of the Davis-Greenstein process should be carried out in the language of Fokker- 
Planck equation (Jones & Spitzer 1967). 

^The word "so-called" is very much in order here, because the mechanical mechanisms, too, provide 
alignment relative to the magnetic line, and their very name, "mechanical" simply reflects the fact 
that these effects are not purely magnetic but involve the grains' mechanical interaction with the 
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(1995 a,b,c,d,e). Lazarian suggested two mechanisms: the cross-over one and the cross-section 
one, and they show themselves in the case of suprathermal grains only. 

The nature of the Gold mechanism is the following. Each collision of the dust particle 
with an atom or a molecule of the streaming gas adds to the particle's angular momentum 
a portion perpendicular to the relative velocity. As explained in one of the above footnotes, 
the major-inertia axis of the body tends to align with the angular momentum. One, hence, 
may say that the interstellar wind will spin-up the granule so that its maximal-inertia axis 
will "prefer" positions perpendicular to the wind. Since the said major inertia axis is, roughly, 
the shortest dimension of the rotator, one may deduce that, statistically, the particles tend to 
rotate with their shortest axes orthogonal to the gas flow. This picture is, though, complified 
by the precession of the magnetic moments (and of the angular momenta that tend, for the 
afore mentioned reason, to align with the magnetic moments) about the magnetic line. This 
mechanism works only for Brownian dust clouds, because it comes into being due to the elastic 
gas-grain collisions to which only thermal granules are sensitive. To be more exact, it is 
assumed here that the precession period is much shorter than a typical time during which the 
grain's angular momentum alters considerably. 

The suprathermally-rotating dust particles ignore the random torques caused by the elastic 
gas-grain collisions, because the timescales for the random torques to alter the spin state 
are several orders of magnitude larger than the average time between subsequent crossovers 
(Purcell 1979). Still, the dust granules do become susceptible to the random torques during 
the brief cross-overs when the granule becomes, for a short time, thermal (i.e., slow spinning). 
This is the essence of the first Lazarian mechanism of alignment, introduced in Lazarian (1995 
d) under the name of "Cross-Over Mechanism". Hence, the first Lazarian mechanism is the 
Gold alignment generalised for suprathermal grains; the generalisation being possible because 
even suprathermal granules become thermal for small time intervals. 

The second Lazarian mechanism, termed by Lazarian (1995 d, e) "Cross-Section Mecha- 



interstellar wind. 
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nism", and studied in Lazarian &; Efroimsky (1996) and Lazarian, Efroimsky & Ozik (1996), 
is not a generalisation of any previously known effect, but is a totally independent, very subtle 
phenomenon. Its essence can well be grasped on the intuitive level: a precessing (about the 
magnetic line) interstellar granule will "prefer" to spend more time in a rotational mode of the 
minimal effective cross section. In other words, the particle has to "find" the preferable mean 
value of its precession-cone's half-angle, value that will minimise the mean cross section. Here, 
the "mean cross section" is the averaged (over rotation, and then over precession) cross section 
of a granule as seen by an observer looking along the direction of interstellar drift. It is crucial 
that, though the alignment is due to gas-grain collisions, it establishes itself not relative to the 
wind direction but relative to the magnetic line about which the spinning grain is precessing. 

Now, the goal is to understand how effective this mechanism is for the dust particles of 
various geometric shapes. Articles (Lazarian and Efroimsky 1996) and (Lazarian, Efroimsky 
& Ozik 1996) addressed the cross-section alignment of oblate and prolate symmetrical grains, 
correspondingly. In the current paper we intend to extend the study to ellipsoidal granules of 
arbitrary ratios between the semiaxes. 

II. STARLIGHT EXTINCTION ON INTERSTELLAR DUST 

A starlight beam passing through a dust nebula gets attenuated. This process is called 
extinction, and it comprises two separate phenomena: scattering and absorption. The final 
result of extinction is a cooperative effect of all dust particles the ray bypasses. 

Concerning the absorption, one may safely take for granted that it is taking place on 
different granules independently from one another. For scattering, though, such independence 
is not, generally, guaranteed. Still, in our study we shall deal with the independent scattering 
solely, in that we shall omit the phase relations between the waves scattered by neighbouring 
grains. This is justified by the starlight not being monochromatic: the lack of coherence in 
it excludes whatever phase-related phenomena^. Thence the intensities of waves scattered 



It may be good to measure the starlight polarisation at separate wavelengths, but such a project 
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from the various granules must be added without regard to phase. Finally, we shall be blithe 
about the multiple scattering, because it brings almost nothing in the observation. Indeed, the 
granules are separated by distances exceeding their size by many orders of magnitude, and the 



To conclude, we shall assume the starlight extinction by different grains to be independent 
and devoid of whatever phase correlations: the appropriate losses in intensity simply add (and 
never return back through radiation transfer). This reasonably simplified approach leaves room 
for a thorough treatment of starlight attenuation by rarified dust. 

III. THE SCATTERING, ABSORPTION, AND EXTINCTION CROSS SECTIONS 

The starlight-scattering cross section on a dust grain is introduced in a pretty standard 
manner. Let us begin, for simplicity, with a monochromatic ray, and then generalise the 
consideration to the natural light. If the incident radiation is a monochromatic plane wave of 
frequency u and intensity Jq, the observer located at point (r, 6*, 0), relative to the scatterer, 
will register intensity proportional to /q and inversly proportional to r~^: 



the wavelength- and angle-dependent factor D having the dimension of squared length. From 
the physical viewpoint, it can have such a dimension if and only if it is proportional to the 
squared wavelength A or, the same, inversly proportional to the squared wave vector k = 
ijjjc = 2n/\ : 



Now let us naively suppose that all the photons entering the scattering region through area 
dCscat get scattered into the solid angle dQ about the direction {9,(j)). If this, purely corpus- 
may require a new, more sensitive, generation of observational means. 
^This is not necessarily true for circumstellar environments where the dust is more abundant, and 
radiation transfer is taking place. 



optical depth of most interstellar^ dust clouds is well below unity. 




(3.1) 




(3.2) 
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cular, interpretation were correct, the energy conservation law would read: 

I{e,(t),uo)dnr^ = lodCscat . (3.3) 

By comparing the latter with the former, we would then arrive to the expression for the 
scattering differential cross section: 

dCscat ^ F{e,<f)) 
The full cross section would then be given by the integral: 



(3.4) 



Cscatiu;) = J dn . (3.5) 

Needless to say, it would be physically incomplete to interpret Cscat simply as the incident 
wavefront area wherefrom the photons get scattered off their initial direction. As well known 
since the times of Newton and Huygens, the corpuscular interpretation neglects the interference 
of the scattered and incident components of light. Therefore, it will, for example, fail to 
correctly describe the forward-scattering. The problem is somewhat subtle. On the one hand, 
the above mathematical expression for differential cross section formally remains correct for 
whatever finite value of the scattering angle 6 . Indeed, for an arbitrarily small, though finite 
angle, the observer potentially can distinguish between the primary and the scattered images. 
To that end, he will have to employ a sufficiently powerful telescope located at a sufficiently 
remote distance from the scatterer. On the other hand, though, the needed resolving power of 
the telescope must be achieved by increasing the size of the object lens. The finite radius of 
lenses thereby imposes a restriction on the scattering angle values for which ( |3.5| ) is meaningful]^. 
Hence, the expression for full cross section Cscat is, in fact, of no physical interest. It 



^Suppose, our telescope is aimed at a distant star, the scatterer being slightly off the line of sight. 
In order for the secondary image to get into the object lens, the scattered photons must be deflected 
at angles not exceeding R/r , with R and r being the radius of the lense and the distance to the 
scatterer. At the same time, the angular resolution of the lense is less than X/R . This results in the 
trivial inequality X/R < R/r , whence R > ^fXr . 
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corresponds to no physical measurement, because it is pointless to carry out the integration 
over too close a vicinity of ^ = . Our is the case: with the distances grossly exceeding 
the device size, the observations of starlight are performed at effectively zero scattering angles; 
so the telescope does not distinguish the forward-scattered light from the primary wave. To 
account for their interference, one has to make a simple estimate. Let the incident radiation 
be described by wave 



Uo 



/o(a;) exp{ — ikz + iuut} 



(3.6) 



Then the scattered wave in the distant zone will read: 

Qm A\ fTT^ exp{-ikr + iut} exp{-ikr + ik z} 
u = S{e,(l)) Jlo{uJ) — = S{e,(j)) — Uo , (3.7) 



ikr ikr 

the relative amplidude S{6,(j)) being, generally, complex. (Evidently, \S{6,(j))\'^ = F{6,(j)) .) 
Consider the case of forward-scattering, when ^ = 0. Let the telescope with object lens' area 
A be located far afield and register a combined image of the incident and forward-scattered 
beams. The aperture plane being z = const, a point of the object lens is located at distance 



r = z + 



2z 



(3. 



from the scatterer. The amplitude at 6' = will be denoted by S{0). The superimposed 
amplitudes will, together, give: 



(3.9) 



the second term being a small correction to the first. Hence, the intensity incident on this 
point of the lens should be: 



\uo + u\ 



ikz 



2z 



(3.10) 



integration whereof over the objective will yield the total intensity of the combined image. As 
well known. 



/ exp (—ikx^/2z^ dx 



I271Z 

ik 



(3.11) 
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wherefrom 



j j \uo + dx dy ~ /o(<^) 



A-^ Re{SiQ)} 



(3.12) 



kz [ i ik 

the integration being carried out over the object-lens area A. Since in the second term we had 
negative-power exponentials, there was nothing wrong in approximating its integral over the 
finite aperture by extending its limits from — oo through oo : for the illustrative estimate, the 
small difference between the Fresnel and Gaussian integrals is irrelevant. 

Evidently, the second, negative term in the latter expression gives the amount by which 
the fight energy entering the telescope is reduced by the scatterer. This reduction is called 
forward-scattering cross-section: 

Cf^M = ^Re{Sm (3.13) 

Roughly speaking, the observer will get an impression that a certain share, Cf^g/-^ of the 
object lens area A is covered up. This shows the fundamental difference between the regular 
scattering cross section and forward-scattering cross section: while the former is associated 
with the area of the incident wavefront, the latter is associated with the area of the observer's 
aperture. This profound difference stems from the fact that the front-scattering fight cannot be 
physically separated from the incident wave. As agreed above, we consider only optically-thin 
clouds. This means that we totally ignore Cgcat , but do take C/_s into account. 

Physically, it is quite obvious that absorption will come into play through adding some 
Cabs to C/_s . Even less light will reach the lense, and the observed intensity will be: 

/o(a;) [A - Cf.sioo) - CaU^)] (3.14) 

In the end of the preceding section we agreed that the light extinction by different granules 
is independent and is free from phase correlations: the intensity losses simply add. This will 
result in the extinction cross sections of the single granules added to give the extinction cross 
section of the entire cloud (for a detailed explanation see van de Hulst 1957, p. 31 - 32). 
Finally, for whatever real observation, the above expression must be multiplied by the window 
function W{ijj) of the device, and integrated over uj . All in all, the resulting attenuation will 
be expressed by the extinction cross-section: 
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(3.15) 




for its cross-sections. If W{uj) carves 



Cext ~ i^f-s + ^Ibs ) 



(3.16) 



All the above is valid for both polarisations independently. So the scattering, forward- 
scattering, absorption, and extinction cross sections may be introduced for the two polarisations 
separately. 

IV. THE MEASURE OF ALIGNMENT, AND ITS RELATION TO OBSERVABLE 



In the Introduction we explained what it means for an interstellar grain to be aligned. For 
the effect to be quantified, it should be endowed with some reasonable measure, one that would 
interconnect the dust dynamics with the starlight-polarisation degree. 

Linear polarisation essentially means that, if the ray propagates in the Zg direction, there 
exist two (orthogonal to Zo and to one another) directions, Xg and i/o , appropriate to the 
maximal and minimal magnitudes of the electric field in this ray. The subscript "o" signifies the 
observer's frame. The question now is: how will these maximal and minimal magnitudes E° 
and Ey (or, equivalently, the maximal and minimal intensities, E°'^ and E"'^ ) evolve along 
the line of sight, within the cloud? Properly speaking, one should talk about the ensemble- 
averages of these intensities: {E°'^) and {Ey"^) , the averaging being implied first over the grain 
orientation (relative to its angular momentum), then over the angular-momentum's precession 
about the magnitic field and, finally, over the half-angle j3 of the precession cone. (The 
direction of magnetic field will be assumed to be constant over the line of sight, within the 
cloud.) While each of the first two averagings will be merely an integration over a full circle, 
the latter averaging will involve some distribution function for f3 . This distribution function 
should be supplied by the detailed theory of a particular orientational mechanism dominating 
the alignment process. 



QUANTITIES. 
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In neglect of the secondary scattering, the decrease in intensity, dl , is proportional to the 
length dz and to the dust-particle density n in the cloud: dl = — Cext n dz , with Cext 
being the afore mentioned extinction cross section ( p.l5| ). What we observe is the intensities 
at the exit from the nebula. Call them {E°'^) and {Ey"^) . Then 

(Ef) ~ exp(-C>/) , (Ef) ~ exp(-C>/) , (4.1) 

/ being the depth of the cloud as seen by the observer, and C° y being the mean extinction 
cross-sections for the two linear polarisations orthogonal to the observer's line of sight. 

Suppose that, prior to entering the nebula, the starlight was unpolarised. One can char- 
acterise the polarising ability of the cloud by the measured flux intensity / as a function 
of the angle of rotation of some analysing element of the telescope. In practice, they rather 
employ a relative measure, Pext , which is the degree of polarisation due to selective extinction 
(Hildebrand 1988): 



p ^ 4 - ly _ {Ef) - {Ef) 
- Ix + Iy~ {Ef) + (Ef) 

As follows from the above formulae, 

exp ( - C° n / ) - exp ( - n / 



(4.2) 



ext 



exp{-C° n I) + exp - C° n I 



tanh 



C°nl - C°nl 



C°nl - C°nl 



(4.3) 



2 

the approximation being valid for Pext ^ 1 (with no need to assume that C° n I <^ 1 ). 
One can also introduce the polarisation per optical depth: 

J^ext _ ^x ^^^^^ 



(Q nl + C°nl)l2 Q + C° 

No matter what measure of polarisation one prefers, this measure involves the difference {C° — 
Cy) that depends on the extinction properties of a single grain and on the degree of alignment 
in the cloud. The topic was addressed by many. A brief conclusion that saves type will be: 
no matter which alignment mechanism dominates, the difference (C° — Cf) should be 
expressed as a function of a single granule's extinction cross sections and of the magnetic 
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field direction (relative to the line of sight). Naturally, the said difference will be also a 
functional of the precession-cone half-angle disribution. This half-angle, often denoted as 
(3 , comprises the angular separation between the magnetic field and the particle's angular 
momentum precessing thereabout. The statistical distribution of (3 over the ensemble depends 
upon the dominating orientational mechanism(s), and its calculation is a technical issue which 
sometimes is extremely laborious and involves numerics. 

The expressions for C° and C° in terms of the afore mentioned arguments were given, 
for the cases of oblate and prolate symmetrical granules, in Greenberg 1968, Purcell & Spitzer 
1971, Lee & Draine 1985, Hildebrand 1988, and Roberge & Lazarian 1999. To fulfil the goal 
of our study, we must generalise those results for the case of triaxial ellipsoid. To this end, we 
introduce the extinction cross sections Cx , Cy , Cz of the grain, for light polarised along 
its minimal {X), middle {Y) and maximal {Z) inertia axesQ. The calculation, presented in 
Appendix A, results in the following relation: 



where 7 is the angle between the magnetic field and the plane of sky (Fig. 1), and R is the 
so-called Rayleigh reduction factor defined as 



f3 being the half-angle of the precession cone described by the grain's angular momentum 
about the magnetic line (Fig. 2). As already mentioned, this angle is not the same for all 



These cross sections characterise particular species of dust. Computation of the granule extinction 
cross sections is comprehensively discussed by van de Hulst 1957. (See also Martin (1974) and Draine 




" Cx + Cy 
. 2 



Cz R cos^ 7 



(4.5) 




(4.6) 



& Lee (1984).) 
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FIGURES 




FIG. 1. The line of sight, the plane of sky, and the direction of magnetic field. Axes y° and 
are chosen to belong to the plane defined by the line of sight and magnetic line. 




FIG. 2. Relative positions of the coordinate systems associated with the magnetic line and with 
the angular momentum. The latter system is obtained from the former through a rotation about axis 
, by angle <p , and a subsequent rotation about axis x'^ , by angle (3 . 
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grains, but obeys some statistical distribution determined by a particular physical setting. The 
above formula for R was obtained in assumption of principal rotation (when the alignment of 
the principal axis with the angular momentum is promptly enforced), valid for suprathermally 
rotating grains (Lazarian & Efroimsky 1999). In the thermal case, when angle 6 between 
the angular momentum and the major-inertia axis is nonzero, the Rayleigh factor would look 
more complicated: 



Returning to (|4.6| ), we would point out the evident fact that, in the absence of alignment, the 
ensemble average of cos^ f3 is equivelent to averaging simply over the solid angle. It entails: 
(cos^/5) = 1/3. This nullifies the R factor and makes the radiation unpolarised: C° = C° . 



V. CALCULATION OF THE RAYLEIGH REDUCTION FACTOR IN THE CASE 
WHEN THE CROSS-SECTION MECHANISM IS DOMINANT 

The basic idea of the cross-section orientational process (pioneered by Lazarian 1995d, 
1995e) comes from the fact that the average time between two sequent crossovers is proportional 
to a typical lifetime of an active site ("Purcell rocket"). Each such site is eventually "poisoned" 
through the evergoing accretion of atoms brought by the interstellar wind. Emergence of new 
active sites leads to cross-overs. Henceforth, the higher the accretion rate, the shorter the 
average lifetime of a typical Purcell rocket. Now, since the atoms adsorbed by the surface get 
delivered through the gas bombardment, one may state that the said lifetime is proportional 
to the rate of gas-grain collisions. The latter rate, in its turn, is proportional to the gas-grain 
effective cross-section, i.e., to the averaged (over the period of precession about the magnetic 
line) cross section of a granule, as seen by an observer looking along the line of gas flow. To 
summarise: the average time between two sequent crossovers is proportional to an active-site's 
lifetime, the latter being proportional to the rate of gas-grain collisions, which in its turn is 
proportional to the effective cross section of the precessing granule in the flow. All in all, the 
dust particle will spend longer times in rotational states of smaller effective cross section. The 
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cross-section mechanism is essential for the both cases of rapid and slow flowsQ. 



This line of reasoning, developed by Lazarian, brings into play several time scales. One 
is the period trot of grain spin about its own rotation axis. The other is the wobble period 
twobbie , i-e., the period of precession of the angular velocity f2 and major- inertia axis Z 
about the angular momentum J . Third is the period T of precession of J about the 
magnetic field B . The fourth one is the mean time r between subsequent flip-overs of the 
granule. 

Suprathermal rotation is swift: time trot is much shorter than the other time scales 
involved. Time t^obUe is irrelevant in the suprathermal case, because in this case we neglect 
the wobble: as explained in Section I, the grain's magnetic moment, the angular-velocity 
vector, and the maxumum- inertia axis are all aligned along the angular momentum, and they 
all precess about B , always remaining parallel to one another. The rate of this precession 
about B is slower than the granule's rotation, but still rapid enough: as mentioned in Section 
I, the precession period T is much shorter than a typical interval r between cross-overs. 
A cross-over of a spinning particle can happen for one (or both) of the two reasons: (1) spin 
damping through collisions with gas atoms, and (2) grain resurfacing that alters positions of 
active cites. Without going into detailed dynamics, let us assume that on the average a cross- 
over takes place after the granule experiences N collisions. Suppose, this amount of collisions 
is achieved during time r : 



^^For high ( > 2 km/s) relative velocities, the cross-section mechanism coerces the grain to align 
in the same direction as the cross-over mechanism does (Lazarian 1995d). For slower flows, the 
stochastic torques produced by the Purcell rockets exceed the torques caused by collisions with gas. 
Thence, no considerable alignemnt will arise during the flip-overs. This makes the role of the cross- 
over mechanism marginal. Therefore, one may expect that the cross-section mechanism will dominate 
in slow flows. Still, the flow should be, at least, mildly supersonic, in order for the drift to dominate 
over the stochastic motion of individual atoms 
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n being the density of atoms, u being the speed of gas flow, and (S)^ being the cross-section 
of the gas-grain interaction, averaged over the grain's (principal) rotation about J : 

(S), = ^ £^ S(^, (3, 0, rj) drj , (5.2) 

the angles $ , /3 , , and rj being as on Figures 2-3. Even after this averaging, (S) 
remains time-dependent, because of the precession of J about B. For times r much longer 
than the precession period T , 



N 



whence 



T = ^ 5.4 



It can be also re-written as 



^- ^ - ^ (55) 

nu (l/27r) /o^- {S)^ d<P nu (l/27r) ( {S),)^ ' ^ ' ^ 

being the angle as on Fig. 2: during the precession period it changes from zero to 27r . The 
usefulness of the latter expression becomes evident if one recalls that the probability to find 
a granule in a certain spin state is proportional to the time it stays there. After averaging 
over the rotation about J , and after a further averaging over the precession, the so averaged 
spin state depends upon two agruments: the precession-cone half-angle (5 , and the angle $ 
between the magnetic line and gas drift (sec Fig. 3). Thence, what the above formula gives us 
is the (not yet normalised) distribution of the (doubly averaged) spin states over (3 (angle 
$ being fixed and playing the role of parameter). What then remains is simply to normalise, 
i.e., to divide r by its integrand over the solid angle. Thus we come to distribution 

with the normalisation constant equal to 

C = r dp sin/3--^ (5.7) 

Jo { {S}r,)^ 
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and the average defined as 

{{S)r,)4> = £^dcf> fj dr^ S (5.8) 
This distribution being on our hands, the Rayleigh reduction factor is straightforward: 




This is where the physics ends and mathematics begins. 

VI. THE DISTRIBUTION OF (3 FOR TRIAXIAL-ELLIPSOID-SHAPED 

GRANULES 

In this section we shall compute the distribution p{(3) which will, in fact, be a function 
of several arguments. One is, naturally, (3 itself. Another will be the angle ^ between the 
magnetic line and the gas flow. Above that, p will depend upon the geometry of grain (in 
assumption that all particles in the cloud are alike). 

To make the section readable, we shall move most part of the calculations to the Appendix. 
Remarkably, these calculations include not just exercises in geometry but also elements of the 
variational calculus, a heavy-duty tool seldom required in astrophysics. 

We shall model the cosmic-dust particle with an ellipsoid of semiaxes a > h > c . Three 
coordinate systems will be employed. Frame (X, Y Z) will be associated with semiaxes 
a, b, c , correspondingly. As we know from Section I, in the suprathermal case the major- 
inertia axis and the angular velocity are aligned with the angular momentum J . So vector 
J is pointing along Z . Another frame, {xb , Vb , zb) will be associated with the magnetic 
field B , axis yB pointing along B (Fig.2). The direction of gas flow will be denoted by 
unit vector f with components {Xf, Yf, Zf) in the body frame. The angle between f and 
axis Z will be called a (Fig. 3). The third coordinate system needed, {xj , yj , zj) 
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FIG. 3. Relative orientation of the gas flow (depicted by vector f ), the 
ellipsoid's principal axes, and its angular momentum J . In the 
principal rotation state the major-inertia axis Z of the body 
(which is its shortest dimension) is aligned with J . 

will be associated with the angular-momentum vector so that yj will be parallel to J and, 
therefore, to Z . Then axes xj and zj will belong to plane (X, Y) . One is free to 
parametrise the rapid rotation of the granule about its major-inertia axis Z by angle r] 
between the least-inertia axis X and the gas flow projection on the plane (X, F) . The 
suprathermal spin about J is a much faster process than the precession of J about the 
magnetic line. Therefore, while averaging over r\ , one may assume the orientation of J 
being unaltered. This means that during several rotations of the granule (about J ) the angle 
a. between J and the gas flow may be assumed unchanged. We shall also need angle $ 
between the gas flow and the magnetic field, and angle /3 between the magnetic field and the 
angular momentum. Finally, we shall parametrise the precession of J about B by angle 
(Fig.2). This angle is constituted by axis zb and axis z' (which is the projection of J on 
the plane perpendicular to B ). Without loss of generality, one can direct axis xb along the 
gas-flow projection on the plane perpendicular to B . 
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As evident from Fig. 4, angles a , (3 , ^ and (f) are not all independent. They obey the 
(proven in Appendix B) relation: 



coscK = cos$ cos/3 + sin$ sin/3 sin0 . 



(6.1) 




,piane orthogonal to 
'the magnetic field 



FIG. 4. Axis is chosen to point along the magnetic hne. Axis is chosen to 
point along the projection of the flow on the plane orthogonal to magnetic 
field B . Auxiliary axis z' points along the projection of the angular mo- 
mentum J on the said plane. 

The gas flow speeding by the eUipsoid defines an elliptic curve bounding area E hatched on 
Fig. 5 and depicted by a thick solid hne on Fig. 6. Evidently, 



^ — \Pmin ^ Pmax\ 



^ Pmin Pmax ; 



(6.2) 



Pmin Pmax being the semiaxes of hatched ellipse. Projection of E on the plane 

orthogonal to the gas flow will give us the cross section S of the granule as seen by the 
observer looking at it along the line of wind. As shown in Appendix B, 



S — TT \ Tlx COS 77 sin a + Uy sin 77 sin a + riz cos a 



(6.3) 



where the body-frame components of the auxiliary vector n = Pmin x Pmax expressed by 
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FIG. 5. 



Gas flow passing by a granule is tangential to its surface in certain points. 
Altogether such points constitute an ellipse whose interior is hatched. Its 

semiaxes are vectors p^m Pmax ■ 




FIG. 6. A granule in the gas flow. The thick solid line denotes cross section S (the 
one that is hatched on Fig. 5). The dashed line is the cross section S of 
the granule relative to the flow. Unit vectors normal to these cross sections, 
n and f , are separated by angle d . 



23 



max Jy 



(6.4) 



and its cyclic transpositions. The components of p^j„ and Pmax depend upon the lengths 
a , b , c of the semiaxes, and upon their orientation relative to the gas flow, i.e., upon angles 
T] and a . Angle a , in its turn, depends upon /? , $ and . All this will, eventually, 
enable us to express 5* via a , b , c , (3 , $ , (p ^ind 1] . After that we shall average over rj 
(i.e., over the principal rotation about J ) and over (p (i.e., over the precession of J about 
B ). It will give us the distribution ( [5. 61 ) over /3 , baring dependence also upon a , b , c and 
$ as parameters. 

Calculation of the components of vector n is presented in Appendix C. Plugging of these 
in (|6.3| ) entails the following expression for the cross section of the grain placed in the flow: 

1 1 



TT 



S = — /i2 I sm" a cos a sm r] cos t] 



+ 
+ 



(52 _ X,) (c2 - Ai) (c2 - A2) (62 - Ai) 
1 1 



(c2 - A2) (a2 - AO 



(a2 - A2) (c2 - Ai) 



(a2 - A2) (62 - Ai) (62 - A2) (a2 - Ai) 



(6.5) 



where 



Pi 



sma COST] 
(a^ — Ai)_ 



sma smi] 



b (62 - A 



-1/2 



cos a 



c (c2 - Ai) 



(6.6) 



P2 



sma cos?7 
a (a? — \o] 



+ 



sma sm?] 
h (62 - A2), 



cos a 



c (c2 - A. 



-1/2 



(6.7) 



Al,2 



Q ± ^Q' - AP R 
2 P 



„ /sina cos?7\2 /sina sinr7\2 /cosa 



(6.9) 



Q 



sma cos?7\2 / 2 2\ /sina sinr^^ 2 



(62 + C2) + 



(c2 + a2) + 



/ COS Oi\'^ ( 2 

\ c 



(a2 + 6^ 



(6.10) 



„ /sina cos?7\2 „ „ /sina sinr7\2 „ „ /cosa\2 „ „ 
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What remains is to average S over rj and (using formula and to plug the inverse 

of ( {S)r,)(f, in the expression (|5.6|) — ( |5.8| ) for the distribution . The latter will give us, 
through ( |5.9| ), the Rayleigh reduction factor as a function of the granule dimensions a, b, c , 
and of the angle $ between the magnetic field and the gas drift. 

This work can be performed only numerically, and must be carried out with a great 
care. The difficulty emerges from the fact that the denominators of the terms in square 
brackets in the above expression ( |6.5| ) for S vanish at certain values of the angles and 
at certain values of a, b, c . Fortunately, this is fully compensated by the multipliers 
/ii /i2 I sin^ a cos a sin r] cos r] \ (which is most natural, for the area of a shadow cast by a 
smooth granule cannot have singularities). Still, when it comes to numerics, the mentioned 
issue demands much attention. 



VII. RESULTS AND THEIR PHYSICAL INTERPRETATION. 

The results of computation are presented on Figures 7-9. As expected, the diagrammes 
on all three pictures show the symmetry that corresponds to the invariance under inversion of 
the magnetic field direction. For another easy check-up, we see that on all the diagrammes the 
Rayleigh reduction factor vanishes in the limit of spherical grains. 

We see that the cross-section Lazarian alignment is intensive for oblate granules (Fig. 7), 
and approaches its maximum in the limit of "fiat fiake" shape. The alignment is maximal 
when the fiow is paraller (or antiparallel) to the magnetic line, or is perpendicular thereto. 
Between these extremes, the alignement goes through zero. To explain this, let us consider the 
simple case of fiate fiake, and recall that the granule, roughly speaking, "wants" to minimise 
its (averaged over rotation and precession) cross-section as "seen" by the fiow. When the fiow 
is parallel (or antiparallel) to the magnetic line about which the grain precesses, the average 
cross section is minimal if the fiake has its precession-cone half-angle 
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■ 0.5 1 1.5 2 2.5 3 3.5 

FIG. 7. The Rayleigh reduction factor R as a function of the angle ^ between 
the magnetic field and the gas flow. The case of oblate symmetrical granules: 
the solid line corresponds to the semiaxes ratio 1 : 1 : 1/10 (flat "flakes"), 
the circle line corresponds to 1:1:1/3, the dashed line corresponds to 1:1:1/2, 
the star line corresponds to 1:1:2/3, the dash-dot line corresponds to 1:1: 0.9 . 

13 close to 7r/2 (and its square cosine close to nil). The R factor will be negative. When the 
drift is orthogonal to the magnetic field, the flake "prefers" to minimise its average cross-section 
by rotating at f3 close to zero (with its squared cosine close to unity). The R factor will be 
positive]^ Therefore, R passes through zero when the angle takes some intermediate value, 
that may be different for different ratios of semiaxes. Contrary to the expectations, though, 
this value, <|)°^''^*<= ^ bares no dependence on the semiaxes' ratio. 



Analytical treatment is possible in the cases of $ = (when radiation-pushed grains follow 
the magnetic line) and of ^ = tt/2 (when the grains are subject to Alfvenic waves or ambipolar 
diffusion). For details see Lazarian & Efroimsky (1996) 
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0.04 




■ 0.5 1 1.5 2 2.5 3 3.5 

FIG. 8. The Rayleigh reduction factor R as a function of the angle ^ between 
the magnetic field and the gas flow. The case of prolate symmetrical granules: 
the solid line corresponds to the semiaxes ratio 1 : 1/10 : 1/10 (elongated "rods"), 
the circle line corresponds to 1:1/3:1/3, the dashed line corresponds to 1:1/2:1/2, 
the star line corresponds to 1 : 2/3 : 2/3 , the dash-dot line corresponds to 1 : 0.9 : 0.9 . 

As evident from Fig. 8, the effect is much weaker in the case of prolate grain^. The 
diagrammes are similar to those of oblate case, and the physical interpretation is the same as 
above. Just as in the case of oblate geometry, the R factor goes through zero at some value of 
$ , that may depend on the semiaxes' ratio. Remarkably, in this case, too, such a dependence 
is absent, and all the curves cross the horizontal axis in the same point <|)P''°'"*'^ . 

Moreover, the angles ^p^°'-°-^^ seems to coincide with (within the limits 

imposed by the calculation error) equals arccos {1 / \f3) . Such miraculous coincidence must 
reflect some physical circumstances that are not evident at the first glance. A straightforward 



^■^This case was studied in (Lazarian, Efroimsky & Ozik 1996). Our Fig. 8 is in full agreement with 
Fig. 3 from that paper. Mind, though, that on Fig. 3 in the said paper there is a slip of the pen: in 
fact, ^ was changing not from through 7r/2 but from 7r/2 through vr . 
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FIG. 9. The Rayleigh reduction factor R as a function of the angle $ between 

the magnetic field and the gas flow. The case of triaxial asymmetrical granules: 
the solid hne corresponds to the semiaxes ratio 1 : 0.9 : 0.2 ; 

the circle line corresponds to 1 : 0.7 : 0.2 ; the dashed line corresponds to 1 : 0.5 : 0.2 ; 
the star line corresponds to 1 : 0.3 : 0.2 ; the dash-dot line corresponds to 1 : 0.2 : 0.2 . 

analytic proof of this "coincidence" , even in the simplest, oblate case, is unavailable. 

The third picture. Fig. 9, accounts for the general case of a triaxial body, never addressed 
in the literature hitherto. Since the triaxial case is somewhat in-between the oblate and prolate 
cases, it is not surprising that the diagrammes have similar form. What is surprising, is that, 
once again, all the curves seem to pass zero in the same point, = arccos {1 j\f3) (and, 
for symmetry reasons, in = tt — arccos {1 /\f3) ). 

Interestingly, the Gold ahgnment of thermal dust fails at the same values of $ (Dolginov 
& Mytrofanov 1976). 

This maddening coincidence makes us suppose that this special value, — 
arccos {1 / \f3) , is completely shape-invariant and is independent from the suprathermal- 
ity degree. So here comes "the arccos (i/V5^)-hypoth esis": No mechanical alignment of 
arbitrarily-shaped (not necessarily ellipsoidal) thermal or suprathermal grains takes place, 
when the magnetic line and the gas drift make angle arccos [1 / \f3) or tt — arccos {1 /VS)- 
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VIII. CONCLUSIONS 



In the article thus far, we have investigated the cross-section mechanism of suparthermal- 
grain aUgnment in a supersonic interstellar gas stream. While the preceding efforts had been 
aimed at the cases of oblate and prolate ellipsoidal grains, in the current paper we studied the 
case of triaxial ellipsoid. We provided a comprehensive semianalytical-seminumerical treatment 
that reveals the dependence of the alignment measure (Rayleigh reduction factor R) upon the 
semiaxes' ratios and upon the angle $ between the magnetic line and gas drift. We provided 
a qualitative physical explanation of some of the obtained results. 

However, the most intriguing result poses a puzzle and still lacks a simple physical explana- 
tion. This is the remarkable shape-independence of the critical value of at which R vanishes 
and the cross-section mechanism fails. For all studied shapes (prolate, oblate, and triaxial with 
various ratios of semiaxes), this critical value is = arccos {1 / \f3) ■ We hypothesise that 
this special nature of the said value of $ is shape-independent. 
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APPENDIX A: RELATIONS BETWEEN THE OBSERVER-FRAME AND 
BODY-FRAME EXTINCTION CROSS SECTIONS 



Our goal is to compute the extinction cross-sections C° and C° for the two hnear 
polarisations orthogonal to one another and to the line of sight, Zo- These, so-called "observer- 
frame" cross-sections should be expressed through extinction cross-sections Cx , C'y , Cz 
appropriate to polarisations along the principal axes X, Y, Z of the granule (with X , Y , 
and Z standing for the minimal-, middle- and maximal- inertia axes, accordingly). 

As an intermediate step, let us first calculate the intensities E° ^ and E" ^ appropriate 
to the linear polarisations, Xg and ijo ■ These observer-frame intensities should be expressed 
through the body-frame intensities Ej^ , Ey , E^ appropriate to polarisations along the 
principal axes. To that end, one has to perform a sequence of coordinate transformations. 

The first step is to express E° and E" through the components E^^y^^ of E in the 
coordinate system associated with the magnetic line. In the observer's frame, axis Zo points 
toward the telescope, while i/o may be chosen to belong to the plane defined by Zg and 
B . Then the magnetic line will lay in the {yoZo) plane. A coordinate system associated with 
the magnetic field (Fig.l) may be defined with i/b pointing along B , and with the x axis 
remaining untouched: xb — Xo- The angle between B and yo (equal to that between zb 
and Zo ) will be called 7 . Hence, 

El = E^ , El = E^ COST + sin 7 . (Al) 

The next transformation is from (x^, , zb) to (xj , , zj) , the latter frame being 
affiliated to the (precessing about B ) angular momentum J of the grain. We choose yj to 
point along J, at the angular separation /3 from (Fig. 2). Vector J describes, about B, 
a precession cone of half-angle (5 ; and so does yj about ys ■ An instantaneous position of the 
rotating frame {xj , yj , zj) , with respect to the inertial one, {xb , ys , Zb) , is determined 
by angle (j) . We see that a transition from (xb , , zb) is two-step: first, we must revolve 
the basis about the yB axis by . This will map axes xb and zb onto axes xj and , 
accordingly. Then we rotate the basis about xj by angle /5 . In the course of a precession 
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cycle of J about B , angle (3 remains unaltered, while (f) describes the full circle. The 
relations between the unit vectors are straightforward: 



z' = COS0 - sin0 (A2) 

z"^ = z' cos (3 — sin /3 — cos cos (3 — sin cos (3 — sin (3 (A3) 

y*^ = y^ cos/5 + z' sin/5 = y-^ cos/5 + z^ cos0 sin/5 - ic'^ sin0 sin/5 (A4) 

x-^ = x-^ COS0 + sin0 . (A5) 

Plugging thereof into the right-hand side of the trivial identity x^ + Ey y^ + Ef — 

E = E^x-^ + E^y-^ + Ei z-^ yields: 

E^ = Ei cos(j) - El sin0 sin/5 - sin0 cos/5 (A6) 

E^ = E^ cos/5 - Eis\n(3 (A7) 

£;f = sine/) + cos(/) sin/5 + E^ cos0 cos/5 (A8) 



Lastly, we take into account the position of the grain itself, relative to coordinates {x'-' , , z"^) . 
As explained in the Introduction (footnote 1), it is legitimate in the suprathermal case to 
assume that the major-inertia axis Z of the particle is aligned with the angular momentum, 
i.e., with axis (for details see Lazarian & Efroimsky 1999). Rotation of the granule is, 
thus, assumed principal and may be parametrised by the angle ip between the least-inertia 
axis X and z'^ (this angle is equal to that between the middle- inertia axis Y and x"^ ). 
This gives: 

El ^ - Ex sinV^ + Ey cosi/j (A9) 
E:^ = Ez (AlO) 
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— Ex cosip + Ey sinip 



(All) 



Combination of all the afore presented transformations results in: 

E° ^ E^ ^ Ei cos (f) - sin (/> sin - E^ sin cos = 

— Ex (— smif] COS0 — costp sin0 cos/3) (A12) 
+ E'y (cosV' COS0 — sin-^ sin0 cos/3) — Ez sin^ sin/3 

and 

El = COS7 + E^ sin7 = 
iJj^^ COS (5 — EI sin /3^ cos 7 + sin (f) + Ey cos sin /3 + cos cos /?) sin 7 



= Ex ( — sin ip sin 7 sin — cos ip cos 7 sin /? + cos ip cos sin 7 cos /?) 

+ i?Y {cosip sin 7 sin0 — sin?/' cos 7 sin/? + sin'?/' cos0 sin 7 cos (3) (^13) 
+ (cos 7 cos (3 + sin 7 cos0 sin/?) , 

whence the averaged-over-the-ensemble intensities read: 

{Ef) -\e],(i + (cos^P)) + \eI(i + (cos^P)) + lEl{sm^p) = 



and 



I {Ej, + E^y + El) + l [1 (El + E^y) - El] I ((cos^/?) - 



^l(El + E'y + El) + (El + E 



Ei 



R 



e: 



o2 



: El sin^7 + - cos^7 (sin^p) + ^ sin'^ j (cos^p)) + 

Ey Q sin^ ^ ^ \ ^ ( ^^^^ ^) ^ \ ^^"^^ ^ ( ^)) 
El (cos'^j (cos^/?) + ^ sin2 7 (sin^/?)^ = 



(A14) 



32 



\ {El + + El) + (\ 



cos^ 7 



\ {Ej, + E'y + El) + (i - 



cos^ 7 



El + El 



El 



R 



(A15) 



Following the established tradition, we single out the so-called Rayleigh reduction factor 

R = z. [ { cos" 13 ) . (A16) 



3 
2 



In ( |A14| ) - ( A15 ) averaging over and implies simply (27r) ^ Jq^ dip Jq^ dcj) , while the 



averaging over /? remains so far unspecified; the appropriate distribution depends upon the 
physics of gas-grain interaction and is calculated in Section V. 

Now, that we have expressed the observable intensities {E°'^) and {E°'^) via those 
appropriate to polarisations along the body axes, we can write down similar expressions inter- 
connecting the extinction cross sections. Since the extinction is merely attenuation of power 
from the incident beam, the contributions to cross-section from the body-axes' directions will 
be proportional to the mean shares of power appropriate to these three axes: 



Cx + a 



Y 



+ 3 



- {Cx + Cy) — Cz 



R 



(A17) 



and 



CI 



wherefrom 



Cx + Cy -\- Cz 



+ 



cos^ 7 



{Cx + Cy 



Cy 



R 



(A18) 



-,{C. 



Cz 



R cos^ 7 



(A19) 



APPENDIX B: ELLIPSOIDAL GRANULE PLACED IN A MAGNETIC FIELD 
AND GAS FLOW. SEVERAL USEFUL FORMULAE. 

Here we prove several formulae used in the text. Let us begin with relations for the angle 
a between the gas-flow direction f and the maximum- inertia axis Z of the dust particle 
(Fig. 4). Dot product of the appropriate unit vectors 
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yj = Yb cos/? + zb sin/3 cos0 + xb sin/3 sin0 (Bl) 

and 

f = Yb cos$ + Zb sin$ (B2) 

leads to formula: 

cos a = cos$ cos/3 + sin$ sin/3 sm(f) . (B3) 
Another important relation is evident from Fig. 3: 

cos a = . (B4) 



On Fig.3, the projection of f on (X, Y) will make an angle rj with axis X , such that 

Xf Yf ^ ^ 

cos?7 sin a = = , sin?] sin a = . = (B5) 

V ^/ + ^/ + ^/ V + ^/ + ^/ 

Formulae (p^ ) and (p5| ) will enable us to calculate the grain's cross section relative to the 
wind. The lines of gas flow, which are tangential to the ellipsoid surface, touch this surface in 
points that altogether constitute a curve. It is the ellipse hatched on Fig. 5. Its area is: 

^ — \Pmin ^ Pmaxl — Pmin Pmax (B6) 

Pmin Pmax being its semiaxes. Projection of S toward the plane perpendicular to the 

gas flow is the cross section S of the granule, as seen by the observer looking at it along the 
line of wind (Fig. 6). In other words, 5* is the "shadow" that the granule casts. Evidently, 

S = Icos^l S (B7) 

9 being the angle between the vector f of the gas flow and vector 

n = Pmin X Pmax (B8) 



orthogonal to the hatched ellipse. As follows from (B4- B5) 



n ■ f \n^ Xf + HyYf + ZA 
cos^ - - 



|n| |f| (S/vr) ,/Xj + Yf + 



TT 

— cos?7 sina + sin?; sina + riz cosa| . (B9) 
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Combining the above with ([B7| ), we arrive to: 

S = n \nx cosrj sin a + Uy sin// sin a + cosa| . (BIO) 

APPENDIX C: CROSS SECTION OF A GRANULE ARBITRARILY ORIENTED IN 

THE GAS STREAM 

Every point (X, Y, Z) belonging to the surface of the elhpsoidal grain obeys equation 

V2 y2 72 

g{X, Y,Z;a,b,c) = — + — + — = l. (CI) 

¥ & 



The points, that constitute the boundary of the hatched figure S on Fig. 5, obey ( |Cll) , along 
with one more relation. That second one is the condition of flow being tangential to the surface 
in these points. Stated alternatively: a normal to the ellipsoid in point (X, y, Z) is given by 
vector (X/a^, 1^/6^, Z/(?) , and the flow must be orthogonal to this normal: 

MX, F, Z- Xj,Yf, Zf) ^ ^ + ^ + ^ = . (C2) 

To find vectors pmin and p^ax pointing from the centre to the closest and the farthest points 
of the ellipse S , one has to employ the variational method: 

^^(?-A9-,A)=0 (C3) 

where l'^ = X^ + + , and A, /i are Lagrange multipliers. The latter equation gives 
us the values of X, Y, Z appropriate to the extremal distances from the origin, in assumption 
that constraints ( Pl|) and ( |C^) hold. Three equations (|C3|), for X, Y and Z , yield: 



2 (a2 - A) ' 2 (62 - A) ' 2 (c^ - A) ^ ' 

for the extremal points where vectors ± Pmin ± p^ax ^nd. To find the values of A , 

plug (PI) in (0): 



x^ y2 2^ 

^ = ^ 2a2 [J- A) + ^ 2&2 (fe/_ A) +^ 2c2 (/-A) = ° " ^^'^ 
This entails: 

35 



where 



Q ± ^Q' - AP R 
2 P 



P 



(C6) 



(C7) 



Q ^ - 



+ & 



Yf 



Substitution of (|C^) in ( |(J1| ) will lead to the expression for /i : 



/x^ = 4 



+ 



(a 



(C9) 



^a(a2-A); VM&'-A); \c{c^-X)) 
Since we have two acceptable values for A , we shall obtain four different values for /i : 



(CIO) 



Xt 



Yf 



Zf 



-1/2 



a (a^ — Ai) 



h (62 - Ai), 



c (c2 - Ai) 



(Cll) 



Xf 



^2 



a (a^ — A9 



+ 



Yf 



-1/2 



(C12) 



b (62 - A2); \c (c2 - As); 
and /is = — /ii , /i4 = — /i2 • Simply from looking at how A and enter expressions ( |C4| ) 
for extremal-point coordinates, we see that a change of sign of fi (with A fixed) corresponds 
merely to a switch from p^^^ (or p„„^, ) to - p„i„ (or - p„„^ , appropriately) . Since 
it is irrelevant, for our purposes, which of the two opposite farmost from the origin points to 
call Pmax ^'^d which to call — p^ax (and, similarly, which of the two closemost to the origin 
points to call p^m which to call — p^ax ) 5 shall take the positive values of /x only. 
As a result, a substitution of Ai and /xi in ( p^ will give us the coordinates of one of the two 
farthest, from the centre, points of the boundary of the hatched ellipse S . Similarly, plugging 
of A2 and fi2 will provide us with the coordinates of one of the two closest points. The chosen 
farthest and closest points will have coordinates {pmax)x,y,z and {pmin)x,y,z , appropriately: 



2 (a2 - Ai) 



max Jy 



2 (62 - Ai) 



) i^Pmax^ 2 



Pi Zf 
2 (c2 - Ai) 



(C13) 
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and 



/i2 Xf 



(^Pmin)y 



P2 Yj 



f 



P2 Z 



f 



(C14) 



2 (a2 - A2) ' 2 (62 - A2) ' '^"""'^ 2 (c2 - A2) 

Further substitution of these expressions in ( |6.4| ) and in its cychc transpositions will lead to: 



Uy Yf = Xj Yj Zj 



Pi P2 




1 






1 




4 


(62 


- A2) (C2 


- Ai) 


(C2 


- A2) (62 


- Ai)J 


Pi P2 




1 






1 




4 


_(C2 


- A2) (a2 


- Ai) 


(a2 


- A2) (C2 


- Ai)J 


Pi P2 




1 






1 




4 


_(a2 


- A2) (62 


- Ai) 


(62 


- A2) (a2 


- Ai)J 



r^zZj = XjYjZj 

Thence the cross section 5* of the grain, relative to the gas-flow direction, will read: 

1 1 



TT 



S = -P1P2 \Xj Yf Zf\ 



+ 

+ 



A2) (c2 - Ai^ 
1 



(c2 - A2) (62 - Ai: 

1 



(c2 - A2) {0? - Ai) (a2 - A2) (c2 - Ai) 
1 1 



(C15) 



(C16) 



(C17) 



(C18) 



(a2 - A2) (62 - Ai) (62 - A2) (a2 - Ai) 

where Ai,2 and ^1^2 are functions of (X/, Yj, Zj) , the latter being functions of angles rj 
and a (while a is, in its turn, depends upon /5 , $ , and cf) ) . All in all, the "shadow" area 
S turns to be a function of angles /? , $ , , and . 
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